Single Mott Transition in Multi-Orbital Hubbard Model 
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The Mott transition in a multi-orbital Hubbard model involving subbands of different widths is 
studied within the dynamical mean field theory. Using the iterated perturbation theory for the 
quantum impurity problem it is shown that at low temperatures inter-orbital Coulomb interactions 
give rise to a single first-order transition rather than a sequence of orbital selective transitions. 
Impurity calculations based on the Quantum Monte Carlo method confirm this qualitative 
behavior. Nevertheless, at finite temperatures, the degree of metallic or insulating behavior of 
the subbands differs greatly. Thus, on the metallic side of the transition, the narrow band can 
exhibit quasi-insulating features, whereas on the insulating side the wide band exhibits pronounced 
bad-metal behavior. This complexity might partly explain contradictory results between several 
previous works. 

PACS numbers: 71. 20. Be, 71.27.+a. 79.60.Bm 
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I. INTRODUCTION 

The nature of the metal insulator transition in multi- 
band materials involving subbands of different widths is 
not yet fully understood. As a result of the geometric 
complexity of many transition metal oxides, the degener- 
acy of valence bands is frequently lifted, giving rise to co- 
existing narrow and wide partially filled bands. A key is- 
sue in these materials is therefore whether these subbands 
exhibit separate Mott transitions or whether single- 
particle hybridization among subbands and inter-orbital 
Coulomb interactions ensure the occurrence of a single 
transition involving all subbands simultaneously. An ex- 
ample is the layer perovskite Sr2Ru04 which consists of 
a wide, two-dimensional d xy band and narrow, nearly 
one-dimensional d xz ,yz bandsi This system is believed to 
exhibit unconventional p-wave superconductivity^ Iso- 
electronic replacement of Sr by Ca leads to a distortion 
of oxygen octahedra and an effective narrowing of the Ru- 
derived t2 g bands, causing a metal insulator transition. 3 
Another case is the layer compound Na x Co02 which as a 
function of doping concentration exhibits a wide range of 
properties, including superconductivity (when hydrated), 
while the parent material C0O2 is an insulator^ 

In previous work^ we investigated the Mott transition 
in Ca2_ 2; Sr 2 ;Ru04 by considering a 3-band model ap- 
propriate for Sr2RuC<4 and varying the on-site Coulomb 
interaction rather than modifying the band width. Ex- 
tending earlier quasi-particle calculations for this system— 
within the dynamical mean field theory^ (DMFT) we 
found that inter-orbital Coulomb interactions lead to 
a significant redistribution of spectral weight between 
t2 g orbitals. As a result, despite the completely dif- 
ferent shapes and widths of the single-particle densities 
of states, for increasing U the d xy and d xz ^ yz subbands 
exhibit similar correlation features, suggesting the ex- 
istence of a single Mott transition. Additional single- 
particle coupling between ti g states caused by octahe- 
dral distortions^ undoubtedly will enhance the trend to- 
wards a single transition. This result was in conflict with 



DMFT calculations by Anisimov et aliS who obtained a 
sequence of Mott transitions as the on-site Coulomb in- 
teraction was increased beyond the single-particle widths 
of the narrow and wide t2 g subbands of S^RuCU, respec- 
tively. 

To avoid the uncertainties stemming from the maxi- 
mum entropy reconstruction-^ of real-frequency spectra 
we recently considered a simple two-band model and fo- 
cused on quantities directly available from the Quantum 
Monte Carlo (QMC) calculations carried out at imag- 
inary times and frequencies^ The subbands were as- 
sumed to have semi-circular density of states of widths 
Wi = 2 eV, W 2 = 4 eV, to be each half-filled, and to 
interact only via local intra- and inter-orbital Coulomb 
energies U and U' — U — 2J, where J is the Hund's 
rule exchange integral. The quantum impurity calcula- 
tions were performed within the DMFT and multi-orbital 
QMC method for T = 125 meV. For increasing U and 
V (with fixed J = 0.2 eV) the subband quasi-particle 
weights Zi were found to vanish at about the same criti- 
cal U c , within a characteristic uncertainty caused by crit- 
ical slowing down. The same behavior was found for the 
quantities Gi((3/2) (fl = l/fc^T) which define the spec- 
tral weight of the coherent subband peaks within a few 
ksT of Ep. Also, two-band calculations at T — within 
the iterated perturbation theory^ (IPT) supported the 
existence of a common transition. These results con- 
firmed our previous finding, namely, that inter-subband 
coupling due to local Coulomb interactions gives rise to a 
single Mott transition rather than a sequence of orbital- 
selective transitions. 

Recently Koga et alii 3 - considered the same two-band 
Hubbard model involving semi-circular densities of states 
with W% — 2 eV, W 2 = 4 eV and on-site Coulomb inter- 
actions specified by U = U' + 2 J. Also assuming both 
subbands to be half-filled, they studied the effect due to 
larger Hund's rule coupling by choosing J = {7/4 and 
U' = U/2. To be able to reach lower temperatures, 
the DMFT calculations were performed within the ex- 
act diagonalization (ED) methocUa and a linearized two- 
site version of the DMFTii^ Surprisingly, both impu- 
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rity treatments indicate the existence of orbital-selective 
metal insulator transitions, i.e., the narrow subband ex- 
hibits a smaller critical U c than the wide subband, im- 
plying a region of coexisting metallic and insulating sub- 
bands in spite of the inter-orbital Coulomb interactions 
defined by U' and J. 

The conflicting results obtained within these various 
approaches underline the subtle and complex nature 
of the Mott transition in multi-orbital systems. Since 
the non-isotropic multi-band character among transition 
metal oxide materials is the rule rather than an excep- 
tion there is evidently a need to clarify this fundamental 
issue. Further theoretical studies are required to test the 
reliability of the DMFT method and various quantum 
impurity treatments. Contradictory results for a rather 
simple model which nevertheless captures a key feature 
present in many real systems obviously shed doubts on 
predictions obtained for more complex materials. 

In the present work we reconsider the same two-band 
Hubbard model as in Refs. 11 and 13. Since the QMC 
calculations are difficult to perform at low temperatures, 
we consider first the simpler iterated perturbation the- 
ory at finite T. Although this scheme is not reliable on 
a quantitative level, it can serve as a useful guide for 
qualitative purposes. At sufficiently low temperatures, 
the IPT reveals a single first-order transition for both 
subbands, confirming our previous results. Because of 
the numerical simplicity of this model we are able to ob- 
tain the entire T — U phase diagram. We then study 
the same model within the two-band QMC method. Al- 
though convergence is considerably more difficult at low 
temperatures, these results also suggest a common Mott 
transition. On the other hand, the interpretation of this 
multi-band phase transition is complicated since in both 
phases the subbands exhibit quite different correlation 
features due to their different single-particle properties. 
Thus, in the metallic phase, the narrow band is much 
more correlated than the wide band and shows a mixture 
of both metallic and insulating features. Similarly, in the 
insulating phase the narrow band exhibits a gap while 
the wide band reveals pronounced bad-metal behavior. 
Thus, the transition is partially incomplete for individ- 
ual subbands. This complexity of the quasi-particle spec- 
tra, and the difficulty of clearly identifying metallic and 
insulating properties at finite temperatures, presumably 
is the origin of some of the contradictions between the 
previous works. 

The outline of this paper is as follows. In the next 
section we discuss the DMFT results derived within the 
IPT as a function of temperature. Section III contains 
the results obtained for the QMC method. Section IV 
provides a summary. 



II. IPT - DMFT 

Let us consider the paramagnetic metal insulator tran- 
sition in a two-band Hubbard model consisting of sub- 



bands of width Wi = 2 eV and W 2 = 4 eV. The 
densities of states are assumed to be semi-circular: 
Pi(to) = 4/(7rWj)(l — 4o; 2 /M / i 2 ) 1 / 2 , corresponding to non- 
hybridizing Bethe lattices. Both bands are taken to be 
half-filled. The interacting Green's functions at imagi- 
nary frequencies are defined as 

Gi(iu n ) = fdu . PM T( . s , (1) 

J lU) n + fl-OJ - hi (tLL) n ) 

where w n >o = (2n + V^nksT are Matsubara frequencies 
and fi is the chemical potential. £j(iw n ) are the subband 
self-energies which must be determined self-consistently. 
Removal of the self-energy from the central site yields the 
impurity Green's functions 

gi(iw n ) = [G^iiun) + Ei(iw„)] _1 . (2) 

The corresponding Fourier transforms at imaginary times 
are denoted as gi{r). Since both bands are assumed to 
be symmetric and half-filled, these functions satisfy the 
conditions <?j(0) = g%{j3) = —0.5 and <7i(r) = gt((3 — 
t). For purely on-site Coulomb interactions specified by 
intra- and inter-orbital matrix elements, the self-energy 
components in second-order perturbation theory are£ 

= f dre^ [<7 lfl ?( T ) + a 2ft (r) 3 2 (r)] , (3) 
Jo 

where j — 1(2) for i = 2(1) and the coefficients at are 
defined as a x = U 2 and a 2 = J 2 + 2{U' 2 + J 2 - U'J). 
Since we take both subbands to be half-filled we do not 
consider here corrections to the bare second-order self- 
energy which should be introduced in the absence of 
particle-hole symmetry^ At half-filling, the first-order 
Hartree-Fock terms are identical for both bands, £f F = 
1.5C7 — 2.5 J, and can be absorbed by the chemical po- 
tential. Throughout this paper we assume J — U/A and 
V = U/2, unless noted otherwise. 

The self-energies Y<i(iLu n ) are now inserted in Eq.(l) 
and the procedure is iterated until self-consistency is 
achieved. The starting self-energy is taken to be zero. 
In particular at low temperatures, the direct and in- 
verse Fourier transforms must be done with sufficient 
accuracy^ Also, at each iteration the symmetry con- 
dition <7i(r) = gi{(3 — r) is enforced to remove spurious 
effects due to numerical inaccuracies. This ensures that 
the Yii(iui n ) remain purely imaginary as they should be 
for symmetric half-filled bands. Typically L = 256 time 
slices and 2 12 Matsubara frequencies were used. Also, 
the real-w integration in Eq. (1) must be performed on a 
rather fine mesh in order to capture the narrow coher- 
ent peak close to the metal insulator transition. This 
can be achieved by linearly interpolating pi(co) and per- 
forming the u> integral analytically between mesh points. 
Once a converged solution is found, the Green's func- 
tions Gi{iuj n ) are used to generate the Fourier transforms 
Gi(r) from which the quasi-particle spectra at real fre- 
quencies are derived via the maximum entropy method, 1 - 
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FIG. 1: Quasi-particle weights Zi of two-band Hubbard 
model as a function of on-site Coulomb energy U for T — 
0.01 eV, calculated within IPT-DMFT. Solid (open) dots: Zi 
(Z-{). (a) increasing U, (b) decreasing U. 
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FIG. 3: Same as Fig. 1 except for (a) T = 0.03 eV and (b) 
T — 0.04 eV. Each Figure shows the hysteresis obtained for 
increasing U (upper points) and decreasing U (lower points). 
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FIG. 2: Same as Fig. 1 except for T = 0.02 eV. 



Figures 1 to 4 show the variation of the subband quasi- 
particle weights Zi « 1/[1 — lmT,i(iu)o) / ujo] as functions 
of U for various temperatures. Up to T = 0.04 eV the 
results show the typical hysteresis behavior obtained also 
for the single-band Hubbard modeliiLi& If U is gradu- 
ally increased from the metallic side the transition occurs 
at a slightly larger critical value than if one begins in the 
insulating phase and gradually decreases U. The hystere- 
sis is most pronounced at low temperatures and becomes 
weaker as T increases. For T > 0.05 eV it is negligibly 
small and the transition begins to resemble the typical 
cross-over behavior. 

The important point of these results is that the quasi- 
particle weights of the narrow and wide subbands ex- 
hibit first-order transitions at precisely the same critical 
U C (T). This picture supports our previous results sug- 
gesting a single Mott transition in a non-isotropic multi- 
orbital environment, in contrast to the orbital selective 
transitions found in Refs. 9 and 11. In line with the nota- 
tion used in the single-band case, we denote the common 
critical value for decreasing U as U c i {T) and for increas- 
ing U as U C 2(T). (These values should not be confused 
with the critical U derived in Ref. 11 for the isotropic 
two-band models.) For T < 0.04 eV the transitions are 
perfectly discontinuous both at U C \{T) and U C 2(T). For 
T > 0.05 eV finite temperature broadening progressively 
dominates the transition region. 

Figure 5 shows the phase diagram obtained from the 
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FIG. 4: Same as Fig. 1 except for (a) T = 0.05 eV and (b) 
T = 0.1 eV. Hysteresis effects are negligible at these temper- 
atures. 
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FIG. 5: Phase diagram for two-band Hubbard model, calcu- 
lated within IPT-DMFT. Open dots: critical Coulomb ener- 
gies U c i(T) and U C 2(T) derived from Figs. 1 to 4. The lines 
are guides to the eye. Their extrapolation yields the common 
second-order critical point marked by a solid dot. The actual 
first-order line U C {T) at low temperatures lies between U c i(T) 
and U C 2(T) and is sketched by the dashed line. 



functions Zi(T,U). The overall shape of the phase 
boundaries agrees well with those obtained for the single- 
band Hubbard model£i£ According to the single-band 
phase diagram, the actual first-order metal insulator 
transition takes place along an intermediate line U C (T) 
such that U c i{T) < U C (T) < U c2 (T) (schematically indi- 
cated by the dashed curve) . Extrapolation of the phase 




FIG. 6: Self-energy T,i(iu3 n ) for two-band Hubbard model 
at T = 0.02 eV near U c2 {T) w 2.91 eV at small Matsubara 
frequencies. Solid (dashed) curves: narrow (wide) subband. 



boundaries suggests a common transition temperature of 
about T c as 0.05 eV. This value agrees well with the 
average between the critical temperatures derived for 
isotropic narrow and wide two-band systems, T c u ~ 
0.033 eV and T c22 w 0.067 eV, respectively. [In the 
single band case the IPT-DMFT yields T c /W » 0.02.1 
According to Eq.(3), for doubly degenerate bands and 
U = 2U' = 4J, this value should be reduced by a factor 

Although inter-orbital Coulomb interactions lead to a 
common metal insulator transition in the narrow and 
wide subbands, it is clear from Figs. 1 to 4 that the 
subband quasi-particle weights differ greatly both in the 
metallic and insulating phases. In this sense, the sub- 
band transitions are partly incomplete: On the metallic 
side, Z\ is consistently smaller than Z 2 because of the 
more pronounced correlations within the narrow band. 
As as result, as will be discussed below, this band can 
exhibit quasi-insulating behavior. Conversely, on the in- 
sulating side of the transition, Z\ drops almost to zero 
while Z 2 drops to a finite intermediate value and gradu- 
ally decreases towards larger U. This 'bad-metal' behav- 
ior of the wide band above the transition is weak at low 
temperatures but becomes more important towards T c . 
There is no indication of a second first-order transition 
in the wide subband at larger U. 

To analyze the subband quasi-particle properties close 
to the transition in more detail we show in Fig. 6 the self- 
energies Ej («<;„) for T = 0.02 eV near U c2 (T). Slightly 
below the transition at U = 2.91 eV the wide band is 
clearly metallic with Y* 2 (iuj n ) ~ lu h at small Matsub- 
ara frequencies. The narrow band exhibits significant 
deviations from this linear behavior. Only at the two 
lowest frequencies Y,i(iuj n ) is roughly linear in to n . At 
larger frequencies it becomes inversely proportional to 
lu u . Thus, this band is in the intermediate region between 
metallic and insulating behavior, exhibiting a narrow co- 
herent peak near Ep nearly separated from pronounced 
Hubbard bands (see below). At U = 2.92 eV, i.e., just 
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FIG. 7: Quasi-particle spectra Ni(w) for two-band Hubbard 
model at T = 0.01 eV near U c i{T). (a) below the transition: 
U = 2.60 eV; (b) above the transition: U = 2.62 eV. Solid 
(dashed) curves: narrow (wide) subband; dash-dotted curves: 
bare densities of states. 



above the transition, the narrow band has become fully 
insulating with Ei(iw„) ~ In contrast, S2(iw n ) 

approaches a constant in the limit of small to n . Thus, 
the wide band still is in the intermediate range between 
metallic and insulating behavior. These results demon- 
strate that in the vicinity of the Mott transition the two 
subbands exhibit a complicated superposition of different 
cross-over behaviors as they switch between the metallic 
and insulating phases. 

The effect of these subband dependent correlations 
on the quasi-particle spectra is illustrated in Fig. 7 for 
T = 0.01 eV just below and above U C \{T). On the metal- 
lic side, the wide band is perfectly metallic. The narrow 
band exhibits a very narrow coherent peak near Ep and 
nearly vanishing weight between this peak and the large 
Hubbard peaks?i2, Thus, the narrow band shows features 
indicative of both metallic and insulating behavior. Simi- 
larly, above the transition the coherent peak in N\ (u>) has 
vanished and a narrow excitation gap is visible. A^(w), 
on the other hand, still has appreciable spectral weight 
near the Fermi energy. According to the phase diagram, 
at sufficiently low temperatures the true first-order line 
U C (T) approaches the boundary U C 2(T). In fact, within 
a consistent dynamical two-site approximation^ U C (T) 
is only slightly smaller than U C 2(T). Fig. 1 shows that 
near this line the 'bad-metal' character of the wide sub- 



FIG. 8: Quasi-particle weights Zi of two-band Hubbard 
model as a function of on-site Coulomb energy U for T — 
0.02 eV and J = 0. Solid (open) dots: Z\ {Z2) for increasing 
U. Symbols (x) in the range 2.2 . . . 2.5 eV: Zi,Z a ~ for 
decreasing U. 



band above the transition is much weaker than close to 
U cl (T). 

The above results suggest that, although there is a sin- 
gle metal insulator transition for both subbands, slightly 
below and above this transition the quasi-particle spectra 
are rather complicated, partially exhibiting 'coexisting' 
metallic and insulating behavior. In fact, in the range 
0.01 eV < T < T c , the transition in the wide band 
might be more appropriately described as metal/bad- 
metal transition. However, this behavior is associated 
with the fact that at finite T there is no clear distinc- 
tion between metal and insulator. Extrapolation of the 
results shown in Figs. 1 to 4 to even lower temperatures 
suggests that the 'bad-metal' behavior in the wide band 
above U C 2{T) eventually disappears and that the gaps 
in both bands open simultaneously at the same critical 
U C (T). 

Before closing this section we briefly address the ef- 
fect of the Hund's rule exchange coupling on the corre- 
lations within the two-band Hubbard model. The rela- 
tively large value assumed so far, J = C//4, gives rise 
to considerable anisotropy and is partly the origin of 
the qualitative differences between the subbands, in par- 
ticular near the transition. The reason for these differ- 
ences can be understood in terms of the self-energy ex- 
pression given in Eq. (3). For J = {7/4: gi/o\ f» 0.5 
while for J = 0: cr^/vi = 2. Near the Fermi level 
the single-particle density of states scales with the band 
width, i.e., g\{(3/2) = 2g 2 {(3/2) before the Coulomb in- 
teraction is switched on. Thus for J = {7/4, t = (3/2 

one finds: £i/E 2 = ( a i9i + oi9\g\) / '(oiffi + Gi9ig\) = 
(8 + 0.5 x 2)/(l + 0.5 x 4) = 3 whereas for J = 0: 
E1/S2 w (8 + 2x2)/(l + 2x4) = 4/3. The subband con- 
tributions to the self-energy therefore differ significantly 
for J = {7/4, but they are quite similar for J = 0, giving 
rise to much more isotropic behavior. 

Fig. 8 shows the quasi-particle weights Z\ as func- 
tions of U for J = at T = 0.02 eV. The compar- 
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ison with Fig. 2 demonstrates that as a result of the 
greater isotropy of the system the correlation induced 
reduction of quasi-particle weight is now quite similar 
for both subbands. In particular, slightly below and 
above the transition both bands show similar metallic 
and insulating behavior, respectively, in contrast to the 
much more complex superposition of quasi-metallic and 
-insulating behavior obtained for J = U /4. As before, 
both bands undergo a first-order metal insulator transi- 
tion at the same critical U . Since there is neither single- 
particle hybridization nor Hund's rule coupling between 
subbands, the common transition is entirely caused by 
the inter-orbital Coulomb interaction U' = U. With- 
out this interaction, we would have Z\(U) » Z 2 {2U) 
(slight differences arise as long as the temperature is held 
fixed). Thus, the critical Coulomb energies of the sub- 
bands would differ by the same factor as the band widths 
Wi. Fig. 8 shows that V enforces a single transition at an 
intermediate U c . As in Fig. 2, the Zi exhibit similar hys- 
teresis behavior: for decreasing U the critical Coulomb 
energy is U C \(T) w 2.19 eV while for increasing U we find 
U c2 (T) w 2.45 eV. 



III. QMC - DMFT 

In our previous workii we considered the same two- 
band Hubbard model as above with W 2 — 2W\ — 
4 eV. DMFT calculations based on the multi-band QMC 
method were performed for a rather high temperature, 
T = 0.125 eV, which corresponds to the cross-over re- 
gion above the critical temperature T c . For J = 0.2 eV, 
the quasi-particle weights Z^ were found to differ signifi- 
cantly at intermediate Coulomb energies due to the dif- 
ferent subband kinetic energy terms. Nevertheless, with 
increasing U the Z% merge again and, within some un- 
certainty associated with critical slowing down, become 
very small at about the same U C (T). 

Fig. 9 shows the quasi-particle weights Zi as func- 
tions of U obtained from analogous QMC calculations for 
J = U/A, U' = U /2. In the cross-over region at high tem- 
peratures (T — 0.125 eV) the larger anisotropy caused 
by the stronger Hund's rule coupling leads to a pro- 
nounced overall broadening of the transition region. This 
applies also to the isotropic two-band models with de- 
generate narrow or wide subbands. Nevertheless, in the 
anisotropic case (12) the Zi approach one another with 
increasing U and exhibit qualitatively similar asymptotic 
behavior. 

The results for T = 0.031 eV shown in Fig. 9(b) ex- 
hibit much sharper transitions in the isotropic as well 
as non-isotropic two-band systems. (In the single-band 
casein T C Q MC /T C IPT w 0.75. Assuming the same ratio 
for the two-band case we estimate from Fig. 5 for the non- 
isotropic model: T ( 9 MC « 0.038 eV.) In the non-isotropic 
case the overall behavior of the quasi-particle weights is 
remarkably similar to the IPT results discussed in the 
previous section. While Z\ drops almost to zero near 
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FIG. 9: Quasi-particle weights Zi of two-band Hubbard 
models as functions of on-site Coulomb energy U calculated 
within the QMC-DMFT for J = (7/4. (a) T = 0.125 eV, (b) 
T = 0.031 eV: Solid dots: isotropic two-band models with 
degenerate narrow bands Wi = 2 eV (11) or wide bands 
Wi — 4 eV (22); open dots: non-isotropic two-band model 
with W 2 = 2Wi = 4 eV (12). Symbols (+) in(b): results 
for T — 0.020 eV. (c) Hysteresis obtained for non-isotropic 
model at T = 0.020 eV: solid (open) dots: narrow (wide) 
band; upper (lower) points: increasing (decreasing) U. Lines 
are guides to the eye. 



U c2 = 2.1 eV Z 2 drops to a finite value and then decreases 
gradually towards larger U. The characteristic kink in Z 2 
becomes even more pronounced at T = 0.02 eV (see + 
symbols). This change in slope of Z 2 (U) near U c2 sug- 
gests that both subbands undergo a common transition 
at the same Coulomb energy and that the wide subband 
exhibits pronounced 'bad-metal' behavior above U c . 
The hysteresis behavior of the Zi for T = 0.02 eV is 
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FIG. 10: Self-energy Ei(ia)„) for two-band Hubbard model 
at T = 0.031 eV near critical (7 at small Matsubara frequen- 
cies, calculated within QMC-DMFT. Solid (dashed) curves: 
narrow (wide) subband. 



shown more detail in Fig. 9(c). While for increasing U 
both Zi exhibit a change of slope near U C 2 ~ 2.1 eV, 
for decreasing U the common transition occurs at U c \ w 
1.8 eV. (The vertical line of dots at 1.8 eV indicates the 
iterations from the lower to the upper branch.) As in 
the case of the IPT, there is no evidence for a second 
transition in the wide subband at larger U. 

As can be seen in Fig. 10, the self-energies J^i(ioj n ) 
close to the critical Coulomb energy derived within the 
QMC method are consistent with those obtained in the 
IPT (see Fig. 6). At U < 2.4 eV the wide band is in 
the cross-over region between metallic and insulating be- 
havior. For U = 2.7 eV T,2(iw n ) becomes inversely pro- 
portional to u) n , i.e., a gap opens up. The narrow band 
undergoes a similar cross-over behavior, except at slightly 
lower values of U. 

Fig. 11 shows the quasi-particle spectra for T = 
0.031 eV in the vicinity of the critical Coulomb energy. 
When the narrow band is about to open a gap, the wide 
band is still fairly metallic. Its own gap would be even 
smaller and is therefore readily filled by finite tempera- 
ture tails extending from the Hubbard peaks towards the 
Fermi level. Only if U is increased to about 2.7 eV is the 
gap in the wide band large enough to not be obliterated 
by this temperature broadening. 

The variation of Zi(U) with increasing U shown in 
Fig. 9(b) is qualitatively similar to the one obtained 
by Koga et al>i£ within the exact diagonalization treat- 
ment of the non-isotropic two-band model: Near U c ~ 
2.6 eV Z\ practically vanishes while Z2{U) remains fi- 
nite but shows a characteristic change in slope. At larger 
Coulomb energies Zi decreases gradually without any 
clear evidence for a subsequent metal insulator transi- 
tion. (U c is somewhat larger than in our QMC results 
since the ED calculations are carried out at near zero 
temperature. According to the phase diagram shown in 
Fig. 5 U C 2(T) increases slightly with decreasing T.) 

Since the ED-DMFT results are basically similar to 
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FIG. 11: Quasi-particle spectra Ni(u>) for two-band Hubbard 
model at T = 0.031 eV, calculated within the QMC-DMFT. 
(a) U = 2.1 eV; (b) U = 2.4 eV; (c) U = 2.7 eV. Solid 
(dashed) curves: narrow (wide) subband; dash-dotted curves: 
bare densities of states. 



our QMC-DMFT calculations, the different conclusions 
reached in Refs. 11 and and 13 seem to stem primarily 
from different interpretations of the underlying physics: 
In our view there is only one first-order transition, but 
with complex quasi-metallic and -insulating subband fea- 
tures just below and above this transition. This picture is 
based on the simultaneous discontinuous changes of the 
subband quasi-particle weights derived within the IPT- 
DMFT and the consistency between the IPT and QMC 
results. On the other hand, if the vanishing of the Zi is 
used as a criterion, the subbands indeed behave very dif- 
ferently. Whereas for the narrow band a critical Coulomb 
energy can readily be identified, Z 2 falls off gradually 
and becomes negligibly small at much higher values of 
U. This metallic 'tail' of Z2 makes it difficult to identify 
a second critical Coulomb energy. According to our QMC 
and IPT results, this metallic tail should diminish at low 
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temperature. It is not clear at present why the ED results 
in Ref. 13 show such a tail although they apply to near 
zero temperature. Also, it is not clear to what extent 
the different subband lie's obtained within the two-site 
DMFT depend on the minimal one-level representation 
of the subband baths. The energy discretization inher- 
ent in both the ED and two-site methods plays a crucial 
role close to the transition when dynamical correlation 
effects on small energy scales are particularly important. 
It would be interesting to perform ED and 2-site DMFT 
calculations in a wider range of temperatures in order to 
analyze the transition region in more detail. 

IV. SUMMARY 

The nature of the Mott transition in multi-band sys- 
tems was investigated within the dynamical mean field 
theory. Although in practice subbands of transition 
metal oxides are usually coupled via single- as well as 
many-electron interactions, we focus here on the effect 
of inter-orbital matrix elements of the Coulomb energy, 
thus neglecting one-electron hybridization. In contrast 
to isotropic systems consisting of identical subbands, the 
metal insulator transition in materials with coexisting 
narrow and wide subbands turns out to be remarkably 
complex. The present results confirm our previous find- 
ing, namely, the existence of a single transition rather 
than a sequence of orbital-selective transitions as the on- 
site Coulomb energy is increased. Nevertheless, at low 
but finite temperatures below the transition the narrow 
band reveals a narrow coherent peak nearly separated 
from the Hubbard bands. Conversely, above the transi- 



tion the wide band shows pronounced bad-metal behav- 
ior. This coexistence of quasi-metallic and -insulating 
spectral features in the vicinity of the transition makes 
the interpretation of the DMFT results non-trivial and 
partly explains the contradictory findings in Refs. 11 and 
13. 

Since it is difficult to perform QMC calculations at very 
low temperatures, we have carried out extensive DMFT 
calculations within the IPT. The phase diagram is con- 
sistent with earlier results obtained for single-band mod- 
els. In particular, the quasi-particle weights show the 
typical hysteresis behavior giving rise to two first-order 
lines U C \(T) and U C 2{T) defining the stability boundaries 
of the insulating and metallic regions, respectively. For 
T < 0.01 eV the results suggest that above the transition 
the metallic states in the wide band should disappear, 
i.e., that the excitation gaps in both bands should open 
at the same U C (T). At intermediate temperatures be- 
low T c the IPT results are consistent with the QMC cal- 
culations: The quasi-particle weights in both subbands 
undergo simultaneous first-order transitions at U c i (T) or 
U C 2(T). These transitions, however, are not complete 
in one or the other of the subbands. For instance, the 
wide band shows appreciable metallicity above the tran- 
sition. Similarly, the narrow band exhibits partly insulat- 
ing Hubbard peaks below the transition. This mixture of 
metallic and insulating behavior diminishes towards low 
temperatures, i.e., purely metallic or insulating phases 
should evolve in both subbands. 
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